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We consider the one-dimensional lattice model of interacting fermions with disorder studied pre- 
viously by Oganesyan and Huse [Phys. Rev. B 75, 155111 (2007)]. To characterize a possible 
many-body localization transition as a function of the disorder strength W , we use an exact renor- 
malization procedure in configuration space that generalizes the Aoki real-space RG procedure for 
Anderson localization one-particle models [H. Aoki, J. Phys. C13, 3369 (1980)]. We focus on the 
statistical properties of the renormalized hopping Vl between two configurations separated by a dis- 
tance L in configuration space (distance being defined as the minimal number of elementary moves 
to go from one configuration to the other). Our numerical results point towards the existence of a 
many-body localization transition at a finite disorder strength Wc- In the localized phase W > Wc, 
the typical renormalized hopping V^^^ = e'"^^ decays exponentially in L as (InV^^'') ~ 

and the localization length diverges as £,ioc{W) ~ {W — Wc)'^^"" with a critical exponent of order 
vioc — 0.45. In the delocalized phase W < Wc, the renormalized hopping remains a finite random 
variable as L ^ oo, and the typical asymptotic value V^'' = e'"^°° presents an essential singularity 
(InVJ^'') ~ —{Wc — W)~'^ with an exponent of order k ~ 1.4. Finally, we show that this analy- 
sis in configuration space is compatible with the localization properties of the simplest two-point 
correlation function in real space. 



I. INTRODUCTION 

Whereas Anderson localization phenomena [l| are rather well understood for a single particle (see the reviews 0-111); 
the case of interacting particles in a random potential has remained much more challenging (see the review Q). In the 
field of disordered fermions, there has been for instance a lot of works on quantum Coulomb glasses (see for instance 
[^-Q and references therein), including the debate on the existence of a metal- insulator transition for interacting 
electrons in two dimensions (see [ol-fTTI and references therein). In the field of disordered bosons, many studies have 
been devoted to the existence and properties of the superfluid-insulator transition (see [T2| - [T7| and references therein) . 
Recently, the idea to reformulate the many-body localization problem as an Anderson localization problem in Fock 
space or in Hilbert space has been very useful |18l - l25| . In particular, this type of analysis has led to the prediction that 
the conductivity of interacting electrons models could exactly vanish in some finite region of parameters in the absence 
of any external continuous bath [2^ [25j . The reason is that conduction mechanisms based on variable-range hopping 
need a continuous bath to locally supply or absorb energy to permit hopping between levels which are not exactly 
degenerate. Since quantum levels are discrete, the many-particle system can fail to be an effective heat-bath for itself. 
Following these ideas, Oganesyan and Huse have proposed that this type of many-body localization transition could 
be realized in some one-dimensional lattice models of interacting fermions [2^ . Unfortunately, the numerical study 
concerning the spectral statistics alone presented in Ref. [261 ] has turned out to be not completely conclusive as a 
result of very strong finite-size effects. Moreover, the analogy with Anderson localization on high dimensional and 
Cayley tree indicates that criteria based on the level repulsion may not be very appropriate in this case [1^. In the 
present paper, we propose to study the existence of a many-body localization transition in the very same model of Ref. 
[26| by studying another type of observable which contains some information on the localization of eigenstates. More 
precisely, we use an exact renormalization procedure in configuration space to compute numerically the renormalized 
hopping between two configurations as a function of their distance in configuration space. We present numerical 
results on the statistical properties of this renormalized hopping that point towards the existence of a many-body 
localization transition at a finite disorder strength. In the localized phase, we measure a power-law divergence of the 
localization length. In the delocalized phase, we measure an essential singularity for the asymptotic renormalized 
hopping. These scaling laws are thus reminiscent of the Anderson localization transition on the Cayley tree, which is 
the simplest example of Anderson transition in a space of infinite dimension, but the values of exponents are different. 

The paper is organized as follows. In Section UTJ we describe how many-body localization models can be studied 
numerically via an exact renormalization procedure ('RG') in configuration space that generalizes Aoki real-space RG 
procedure for Anderson localization one-particle models. In Section IIIIl we present our numerical results for a one- 
dimensional lattice model of interacting fermions. In section HVl we discuss the similarities and differences with the 
scaling laws of Anderson localization on the Cayley tree, and propose a specific form of finite-size scaling analysis. In 
section |Vl we present independent numerical results concerning the simplest real-space two-point correlation function, 
to test the compatibility with our results obtained in configuration space. Our conclusions are summarized in section 
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II. EXACT RENORMALIZATION PROCEDURE IN CONFIGURATION SPACE 

In this section, we briefly summarize the Aoki real space renormalization ('RG') method for one-particle localization 
models before we describe its generalization for many-body problems. 

A. Reminder on Aoki real space RG for one-particle localization models 

For Anderson localization models, there exists an exact re al-sp ace renormalization procedure at fixed energy E 
which preserves the Green functions of the remaining sites |27H33 |. 

The renormalization (RG) procedure can be applied to any Anderson localization model of the generic form 

i7 = ^l/,,,K><i| (1) 

where Vi,i is the on-site energy of site i, and where Vij is the hopping between the sites i and j. Upon the elimination 
of site io in the Schrodinger equation at energy E 

E^ito) = K:o,.o^(«o) + E ^^o,Mj) (2) 
the remaining sites satisfy the Schrodinger equation at energy E with the renormalized parameters 

-trnew T/ I ^i,ia J /o\ 

- ^i,3 + p _ T. I'iJ 

These renormalization equations are exact since they are based on elimination of the variable V'(*o) in th^ Schrodinger 
Equation. As stressed by Aoki [27I, ^2§\, the RG rules preserve the Green function for the remaining sites. This 
means for instance that if external leads are attached to all surviving sites, the scattering properties will be exactly 
determined using the renormalized parameters (see (32l | for more details). In particular, the renormalized hopping 
between the last two surviving sites (after all other sites have been decimated) determines the two-point Landauer 
transmission between leads attached to these two points [s^ : it decays exponentially with the distance in the localized 
phase, it remains finite in the delocalized phase, and at the critical point it becomes multifractal. 

B. Generalization in configuration space for many-body localization models 

The above RG procedure has the following natural generalization for many-body models. Let us denote C a 
configuration of the many-body problem to write the Hamiltonian as 

Hmany = ^ ^C. ,C, |C^ >< | (4) 

Then the Schrodinger equation projected onto the configuration Cig 

allows to eliminate ^l){Cig). The remaining configurations satisfy the Schrodinger equation with renormalized param- 
eters 

Vr- C ^C- C 

VPT =Vc c + ° ° ' (6) 

These rules in configuration space have been already used for the two-particle ID Anderson tight-binding model [33|. 
In the following, we apply them to a model of interacting fermions that we now describe. 
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C. Application to the interacting fermions model of Ref. f26|| 

In numerical studies of quantum problems containing both interactions and disorder, it is natural to consider first 
the spatial dimension d = 1. The simplest model is then a chain of spinless fermions with nearest- neighbor interaction 
and on-site disorder (see for instance [33 - [36| and references therein) . This type of model can also be studied in the 
language of quantum spin chains as in Ref. [STj . where a powerful time-dependent Density Matrix Renormalization 
Group method has been used to characterize the many-body localized phase. In this paper, we consider the same 
class of model, but with second-neighbor hopping in addition, as in Ref. ^] . More precisely, the model of Ref. [l^l 
is defined by the following Hamiltonian on a one-dimensional lattice of L sites with periodic boundary conditions 



(7) 



with the usual notations : 

(i) Tii — c\ci represents the number of spinless fermion on site i and can take only two values (0 if the site is empty 
or 1 if the site is occupied). The many-body Hilbert space has thus for dimension 

Ml - 2^ (8) 
The spinless character has been chosen to reach bigger sizes L for a given value of the Hilbert space dimension Ml 

(ii) the on-site energies Wi are independent Gaussian variables with zero mean and variance , i.e. W measures 
the disorder strength (we have not used the 'microcanonical constraint' of Ref. [1^ consisting in the requirement that 
(1/L)^w^ should be exactly W"^ to reduce statistical uncertainties). 

(iii) the nearest-neighbor interaction is chosen to be ^ = 2, the hopping terms between nearest-neighbors and 
second- neighbors are chosen to be t = t' = 1. The second-neighbor hopping is included to have non-integrability at 
zero randomness, see [38| for more details on the properties of the model in the zero-disorder limit. 

(iv) the total number of particles is conserved : we study the case of half-filling with L/2 particles for L sites as in 

[2q . The dimension of the Hilbert space is then given by the binomial coefficient 

Physically, the important point is that at leading order, it st ill g rows exponentially in L. 

In summary, we consider in this paper the model of Ref. [26| with exactly the same values of parameters, but we 
study another observable to detect the possible many-body localization transition. We have applied the RG procedure 
in configuration space described above, to obtain, in each disordered sample of even size L, the renormalized hopping 
Vl at some energy E between two configurations Ca, and Cb after the decimation of all other configurations 

Vl^Vc,,cAE) (10) 

We have made the following choices : 

(cl) we consider the zero-energy case E = 0, because it represents the center of the many-body energy levels. 
Indeed, in Anderson localization models, it is well known that energy levels near the center are the more favourable 
to delocalization : if these states are localized, one expects that all other states will also be localized. 

(c2) we have chosen the following configurations : the configuration Ca has all even sites occupied and all odd 
sites empty, whereas the configuration Cb has all odd sites occupied and all even sites empty. Their distance in 
configuration space is thus L/2 (the minimal path to go from configuration Ca to configuration Cb requires L/2 
elementary moves). In the absence of disorder, the model is known to be conducting (see [S^ for a detailed study of 
conductivity properties) ; the two configurations Ca and Cb are equivalent up to a translation of one lattice site, and are 
thus expected to be connected by a finite renormalized hopping. In the presence of disorder, these two configurations 
are not equivalent anymore, and one expects that the renormalized hopping Vl will become exponentially small in L 
for sufficiently strong disorder. 

This choice (c2) of alternate configurations Ca and Cb can be questioned in various ways. For instance, if one 
wishes to maximize the distance in Fock space, one obtains the configurations where all particles are on the first half 
or on the second half : physically, it is however clear that these two configurations are not typical because they are 
extremely inhomogeneous and because only 4 particles (2 particles at each boundary of the macroscopic cluster) can 
move for arbitrary large L (instead of an extensive number of particles for typical configurations) . More generally, in 
contrast to usual Anderson localization models where all sites are equivalent, a new difficulty that arises in many-body 
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localization models is that all configurations are not equivalent : configurations have different hopping connectivities, 
and different interaction energies, so that the configuration space has already an inhomogeneous structure even before 
the introduction of disorder variables. Since an extensive study of the renormalized hoppings in this complicated 
configuration space is not really possible, we have decided to consider only the choice (c2) of alternate configurations 
Ca and Cb in the remaining of this paper. However, to show that our results are meaningful and do not really depend 
of our precise choice (c2), we present in section fVl independent numerical results concerning the simplest real-space 
two-point correlation function. 



III. 



NUMERICAL RESULTS FOR THE INTERACTING FERMIONS MODEL OF REF. n 



In this section, we describe our numerical results concerning the statistics oi Vl B,t E = between the alternate 
configurations Ca and Cb described above for even sizes 4 < L < 12 with corresponding statistics 15.10'' < ng{L) < 
3650 of disordered samples (we have also data corresponding to L — 14 with Ug = 100 samples, but this statistics has 
turned out to be insufficient for most purposes). 



Analysis of the localized phase 
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FIG. 1: (Color on line) Statistics of the renormalized hopping Vl in the localized phase (data for the disorder strength W = 20) 
: (a) Probabihty distribution PL(lnV'i,) of the logarithm of renormalized hopping Vl for the sizes L — 4,6,8, 10, 12. (b) same 
data for the rescaled variable x of Eq. [Tl]for the sizes L = 6, 8, 10, 12 : the convergence towards a fixed rescaled distribution 
P{x) is rapid (we have only excluded the smallest size L = 4 that was too different). 



For strong disorder, we find that the renormalized hopping Vl introduced in Eq. [TU] flows towards smaller and 
smaller values as L increases. As an example for W — 20, we show on Fig. [T]the probability distributions PL{lnVL) of 
the variable In over the disordered samples of a given size L = 4,6,8,10,12 : the regular shift of these histograms 
towards smaller values is clear. We show on Fig. [T] (b) the same data for the rescaled variable 



IuVl-IuVl 



(11) 



where In Vl is the averaged value and where A^, is the width of the distribution PL{\nVL)- One can see on Fig. [l](b) 
that the histograms of the rescaled variable x coincide within statistical fluctuations for L = 6, 8, 10, 12 (we have only 
excluded the smallest size L = 4 that was too different) : this shows that the convergence towards a stable rescaled 
distribution P{x) is rapid for this model. 

We show on Fig. [2] (a) the decay with L of the disorder-average In Vl for various disorder strengths in the range 

6 <W < 20 : these curves correspond to an exponential decay of the typical value V^^^ = e''^ with respect to the 
distance {L/2) in conflguration space 



in(v^r) 



\nVL{W>W,) 



{L/2) 



(12) 



5 




FIG. 2: (Color on line) Exponential decay of the typical renormalized hopping Vj^^^ = e'"^^ in the localized phase : (a) Linear 
decay of In Vl as a function of L (see Eq. I12|) . (b) Behavior of the slope l/^ioc{W) (inverse of the localization length ^ioc{W) 
of Eq. I12|l as a function of the disorder strength W 



where ^ioc{W) represents the localization length that diverges at the delocalization transition 

6oc(W^) ^ {cst){W -Wc)-"'"^ (13) 

We show on Fig. [2] (b) our numerical result for the slope l/6oc(W^) as a function of the disorder strength W in the 
region 6 < W < 20 (below W = 6 we cannot estimate the linear slope anymore). A three-parameter fit of the form 
of Eq. [13] yields a critical point in the range 

5.2 <Wc< 5.9 (14) 

and a critical exponent around 

t^ioc ^ 0.45 (15) 



B. Analysis of the delocalized phase 



For weak disorder, we find that the renormalized hopping Vl introduced in Eq. [TU] remains a finite random variable 
Voo finite as L increases. As an example for W = 2, we show on Fig. [3] the probability distributions Pl^IuVl) of the 
variable In Vl over the disordered samples of a given size L = 8, 10, 12 (for clarity we have excluded the smallest sizes 
L = 4, 6 that were too different) : it is clear that these histograms coincide up to statistical fiuctuations. This should 
be compared with Fig. [T] (a) corresponding to the localized phase for W = 20. In the delocalized phase, the typical 
renormalized hopping = e'"^°° thus remains finite 

\nVL{W < Wc,L) ~ \nVoo{W < Wc) finite (16) 

We show on Fig. [3](b) our numerical estimates of the asymptotic value In Voo(W^ < Wc) as a function of W. We find 
that our data are compatible with an essential singularity behavior of the typical asymptotic hopping 

\nV^P{W < Wc) = InV^iW < Wc) ~ -{cst){Wc - W)-^ (17) 

A three-parameter fit of this form yields a critical point in the range 

5.5 <Wc< 5.7 (18) 

and an essential singularity exponent around 

K~1.4 (19) 

Essential singularities in transport properties have already been found in various disordered models, in particular 
in Anderson localization on the Cayley tree (see (4TI . [45j and references therein) and in superfluid-insulator transitions 
of disordered bosons (see [l^ - [l^ . [Tq and references therein). 
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FIG. 3: (Color on line) Statistics of the renormalized hopping Vl in the delocalized phase (a) For the disorder strength W — 2, 
the probability distribution Pz,(ln Vl) of the logarithm of renormalized hopping Vl remains the same for the sizes L = 8, 10, 12 
(we have excluded the smallest sizes L — 4,6 that were a bit different). This should be compared with Fig. [T](a) corresponding 
to the localized phase, (b) Behavior of the typical asymptotic renormalized hopping V^^ = e'"^°° : InV^^ — InVoo as a 
function of the disorder strength W . 



C. Conclusion of the numerical study 




FIG. 4: Critical behaviors obtained for a critical point at the value Wc = 5.6 (a) Localized phase W > Wc : the plot of 
ln(l/^ioc(VK)) as a function of ln(M^ — Wc) corresponds to the slope vioc — 0.45 (see Eq. I13p (b) Delocalized phase W < Wc : 
the plot of ln(— In Voo) as a function of ln(Wc — W) corresponds to the slope k ~ 1.4 (see Eq. I17|) 

In summary, our numerical data are compatible with a many-body localization transition for the model of Eq. [T] 
The global best value for the critical point seems to be 

Wc ~ 5.6 (20) 

which is somewhat smaller than the critical value suggested by the level statistics study of Ref. [1^. A possible 
reason for this slight difference could be that the level statistics study of Ref. [l^ is based on all levels of all energies, 
that could mix contributions of various types of states (delocalized, localized, and critical), whereas we have chosen 
to work at the fixed energy E = (center of the many-body energy levels). Anyway, taking into account the large 
uncertainties on Wc as estimated from small system sizes, we feel that the two studies point towards the same region 
of disorder strength W. 

For the value of Eq. [20l we show the log-log plots of the critical behaviors on Fig. ID On Fig. |4] (a), we show 
the divergence of the localization length ^ioc{W) in the localized phase W > Wc : the slope vioc — 0.45 (see Eq. 
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[T3)) . On Fig. |4] (b), we show the essential singularity of the typical asymptotic hopping in the delocalized phase 
W < Wc : the slope corresponds to the exponent k ~ 1.4 (see Eq. [T71). Of course, these values are not expected to 
be precise, since they have been obtained from small system sizes and some fitting/extrapolation procedures from the 
raw data. Nevertheless, the emergence of reasonable scaling laws is encouraging. In the following section, we discuss 
the similarity with the scaling laws that appear for Anderson localization on the Cayley tree. 



IV. DISCUSSION : SIMILARITIES AND DIFFERENCES WITH ANDERSON LOCALIZATION ON 

THE CAYLEY TREE 



A. Analogy with Anderson localization on the Cayley tree 



As recalled in the introduction, the reformulation of the many-body localization problem as an Anderson localization 
problem in Fock space or in Hilbert space has been very useful [18l - t25j . The idea is to analyse whether there exists an 
Anderson localization in configuration space, and to study the consequences for real-space properties. The geometry 
of configuration space is usually very different from the regular finite-dimensional lattices considered in Anderson 
one-particle localization models, and has been argued to be qualitatively similar to the Cayley tree [l8l ll9l [2ll . [23 |. 
Since Anderson localization on the Cayley tree has been studied for a long time as a mean-field limit 39^-'45|, results 
and methods have been then borrowed to analyse many-body localization properties of quantum dots JJi, lU , [2l|, . 
This approximation by a tree structure has been however sometimes criticized [l^l . Indeed, in many-body localization 
models, the Fock space or Hilbert space is never exactly a tree, and thus the approximation by a Cayley tree has been 
proposed as a simplifying approximation to obtain an exactly solved model [18J (note however that in [24| , it has been 
argued that an effective Cayley tree structure should actually well capture the properties of low-dimensional electronic 
models). But independently of the technical convenience of the tree structure, we believe that the physically important 
property in this analogy is the 'infinite-dimension' property, defined as the exponential growth of the configuration 
space A/£ with the real-space linear size L 

Ml oc e^"''^^ (21) 

(whereas in finite dimension d, the configuration space of a single particle grows as a power- law Ml oc L"^). As argued 
in [4^ , it is the exponential growth of Eq. [5T] which is directly responsible for the presence of essential singularities 
of transport properties, whereas finite-dimensional lattices are characterized by power-law singularities. From the 
point of view of Anderson one-particle models, the Cayley tree is thus rather 'pathological' since d = 00 turns out to 
be a singular point, and the upper critical dimension is considered to be dc — +00 [46|. From the point of view of 
many-body localization however, the exponential growth of Eq. [5T]is the rule (see for instance Eq. [5]), and thus the 
scaling behaviors that appear on the Caylee tree are instructive, as an example of Anderson localization on a space 
of infinite dimension. In particular, this analysis suggests some specific form of finite-size scaling as we now recall. 



B. Specific form of finite-size scaling in the critical region 



As in our recent study of the Landauer transmission for Anderson localization on the Cayley tree [45| , it is natural 
to assume some finite-size scaling in the critical region of the form 

In (T4^) --L^G (^L^/'^'^Wc - W)) (22) 

where the finite-size scaling exponent vps is different from the localization length exponent vioc- (This is in contrast 
with the scaling theory of localization in finite dimension d, where the finite-size scaling is governed by t'/oc)- 
The matching of Eq. [52] with the localized phase (see Eq. [T^] and Eq. [T5|) yields 

Vloc = (1 - P^FS (23) 
and the matching with the delocalized phase (Eq. [T71) yields 

K = pi^FS (24) 
By consistence, the finite-size correlation length exponent vps is then given by 



VPS = ^loc + K 



(25) 
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FIG. 5: (Color on line) Finite-size scaling analysis of our numerical data for the sizes L = 8 (triangles), L = 10 (squares) and 
L — 12 (circles) according to the form of Eq. [22] with the values p = 0.76 (see Eq. I28p and ufs = 1-85 (see Eq. I25p : this 
plot is rather convincing at criticality and in the localized phase W > Wc, whereas stronger corrections to scaling seem to be 
present in the delocalized phase W < Wc- 



In an exactly solved travelling/non-travelling phase transition where the same type of finite-size scaling occurs [471 . 
the physical interpretation of the finite-size scaling exponent vps is that it governs the relaxation rate towards the 
finite value in the non-travelling phase. For Anderson localization on the Cayley tree, we have checked that this 
interpretation holds (45j |. For the present many-body localization transition, this property cannot be checked with 
our numerical data limited to small sizes. 

Exactly at criticality, we thus expect the following stretched exponential decay of the typical renormalized hopping 



\nVLiWc)~-LP (26) 
where the exponent p is related to the other exponents by (see the scaling relations of Eqs [23] and [24] ) 

— (27) 

From our previous estimates of the exponents vioc — 0.45 (Eq. [T5)) and k ~ 1.4 (Eq. [T9|) . this would correspond to a 
numerical value of order 



p ~ 0.76 (28) 

We show on Fig. [S] the finite-size scaling analysis of our numerical data according to the form of Eq. [52] with the 
values p = 0.76 and vps = 1-85 obtained by consistency from our previous estimates of vio^ and k : the data collapse 
seems satisfactory at criticality and in the localized phase W > Wc, whereas stronger corrections to scaling seem to 
be present in the delocalized phase W < Wc- 

In summary of this discussion, we propose that the scaling laws of many-body localization transitions are generically 
similar to the scaling laws observed for Anderson localization on the Cayley tree, as a consequence of the infinite- 
dimension property of Eq. 1211 However, besides this qualitative analogy, one should not expect an exact equivalence 
in general, and in particular the critical exponents {viocVps, k, p) are not expected to be the same as those of the 
Cayley tree. 



V. NUMERICAL RESULTS CONCERNING THE SIMPLEST REAL-SPACE TWO-POINT 

CORRELATION FUNCTION 

As recalled in the introduction, the idea that many-body localization actually occurs in configuration space is 
very useful, and in this paper, we have adopted this point of view : we have focused on the renormalized hopping 
between two configurations separated by a given distance in Fock space, with the hope that the signatures of the 
transition would be clearer for this observable. Nevertheles, it is of course very important to understand what are 
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FIG. 6: (Color on line) Exponential decay of the typical real space two-point correlation function C*^^{L) = e'"*-^'^' (see Eq. 
1291 ) in the localized phase : (a) Linear decay of InVl as a function of L (see Eq. I12|l . (b) Behavior of the slope l/^c{W) 
(inverse of the localization length ^c{W) of Eq. I30p as a function of the disorder strength W (circles), as compared to l/£,ioc{W) 
(squares) found previously for the renormalized hopping in configuration space (see Eq. [12] and Fig. [2]). 



the consequences of this localization occuring in configuration space for real-space properties. In this section we thus 
present direct calculations of the simplest two-point correlation function 

C{L) = I < l^rmd\c\cL\i^„M > |' (29) 

where \4'mid > is the eigenstate (obtained via exact diagonalization) of the Hamiltonian of Eq. [7] with free ends at 
i — 1 and i — L (no periodic boundary conditions, so that these two points are at distance L in real space), whose 
eigenvalue Emid is in the middle of the many-body energy levels (this energy E„iid fluctuates from sample to sample 
but remains close to the central value E = chosen in (cl) of section Hi Cp . 

We show on Fig. in](a) the decay with the distance L of the disorder-average In C(L) for various disorder strengths 
in the range 6 < < 20 : these curves correspond to an exponential decay with L of the typical value C^^^ = e'"*-^'^) 

HCT)^^^Cl{W>W,) ^ --^ (30) 

where represents the correlation length that diverges at the delocalization transition 

ic{W) ~ {cst){W -W,)-"^ (31) 

We show on Fig. |6] (b) our numerical result for the slope 1 /^c (W) as a function of the disorder strength W in the 
region Q <W < 20, as compared to \/£,ioc{W) found previously for the renormalized hopping in configuration space 
(see Eq. [12] and Fig. [2]). Our conclusion is that up to a numerical prefactor, these two correlation lengths seen either 
in the renormalized hoppings in configuration space, or in the two-point correlation function in real space, contain 
essentially the same information. In particular, a three-parameter fit of the form of Eq. I31l vields values for the critical 
point Wc and for the critical exponent uc ^ ^loc that are compatible with the values estimated previously from the 
data in configuration space. 

In the delocalized phase however, it is not clear to us what are the theoretical expectations for the decay in L 
of the two-point correlation function of Eq. 1291 and our numerical data are not sufficiently clear by themselves to 
indicate which procedure should be used to fit the data in order to obtain information on the critical behavior in the 
delocalized phase. Further work is needed to clarify this point, or to find other real-space observables that display a 
clearer behavior in the delocalized phase. 

To summarize this section, our numerical data concerning the real-space two-point correlation function of Eq. 
[29] indicate that the correlation length measured previously in configuration space is essentially equivalent to the 
correlation length measured in real space. In particular, this shows that the results obtained in configuration space 
do not depend to much on the particular choice (c2) of alternate configurations made in section [IlCl 
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VI. CONCLUSION 

In this paper, we have proposed to study many-body locaUzation transition via an exact renormaUzation procedure 
in configuration space that generalizes the Aoki real-space RG procedure for Anderson localization one-particle models. 
For the one-dimensional lattice model of interacting fermions with disorder studied previously by Oganesyan and Huse 
[26j , we have studied numerically the statistical properties of the renormalized hopping Vl between two configurations 
separated by a distance L in configuration space. Our numerical results are compatible with the existence of a many- 
body localization transition at a finite disorder strength of order Wc ~ 5.6. In the localized phase W > Wc, we have 
found that the typical renormalized hopping V^^^ = e'"^^ decays exponentially in L as (In Vl^'') ~ ~^7^ ^^^d that the 
localization length diverges as ^loc ^ (W — Wc)'"'"" with the critical exponent of order vioc — 0.45. In the delocalized 
phase W < Wc, we have found that the renormalized hopping Vl remains a finite random variable Voo as L — >■ oo, 
and that the typical asymptotic value V^p = e'"^=° presents an essential singularity {InV^''') ^ —{Wc — W)^'^ with 
an exponent of order k ^ 1.4. We have argued that the analogy with Anderson localization on the Cayley tree 
is important as an example of Anderson transition on a space of infinite dimension (in the sense of Eq. I2ip that 
presents essential singularities and that it suggests a specific form of finite size scaling that we have tested. Even if 
the numerical values of the exponents are not expected to be precise, as a consequence of the limited system sizes 
studied L < 14, we hope that the scaling laws that emerge are valid. Of course, it would be very useful in the future 
to test these results with other numerical methods like the Density Matrix Renormalization Group method, that allow 
to study these interacting one-dimensional models for much bigger system sizes [33 - [37l | . Finally, we have shown that 
the present analysis in configuration space is compatible with the localization properties displayed by the simplest 
two-point correlation function in real space. 

In conclusion, the reformulation of many-body localization problems as Anderson localization models in configu- 
ration space raises the question of Anderson localization on specific networks (see [48| and references therein), that 
are completely different from the regular lattices that have been considered in the field of one-particle models. For 
a many-body problem defined on a domain of size L"^, the number of configurations (i.e. the nodes of the network) 
grows exponentially A/l oc e'-'^"*^^ . Each configuration has a different connectivity in this space of configurations, but 
it is typically of order L'^ (assuming a finite density of fermions, with a finite number of short-range hopping for each 
fermion). Besides the interest in specific many-body models, an important issue is of course to understand which 
properties of this complex network are relevant to determine the universality class of the corresponding Anderson 
transition. 
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